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Abstract. The luminosity evolution of stars with highly 
condensed cores surrounded by nuclear-burning shell(s) is 
analytically investigated with the aid of homology rela- 
tions. With respect to earlier works using a similar ap- 
proach (e.g. Refsdal & Weigert 1970; Kippenhahn 1981), 
the major improvement is that we derive all the basic de- 
pendences (i.e. on core mass, core radius, and chemical 
composition) in a completely generalised fashion, then ac- 
counting for a large range of possible physical properties 
characterising the burning shell(s). Parameterised formu- 
las for the luminosity are given as a function of the (i) 
relative contribution of the gas to the total pressure (gas 
plus radiation), (ii) opacity source, and (iii) dominant nu- 
clear reaction rates. 

In this way, the same formalism can be applied to 
shell-burning stars of various metallicities and in differ- 
ent evolutionary phases. In particular, we present some 
applications concerning the luminosity evolution of RGB 
and AGB stars with different chemical compositions, in- 
cluding the case of initial zero metallicity. It turns out 
that homology predictions provide a good approximation 
to the results of stellar model calculations. 

Therefore, the proposed formalism is useful to un- 
derstand the possible differences in the luminosity evo- 
lution of shell-burning stars within a unified interpreta- 
tive framework, and can be as well adopted to improve 
the analytical description of stellar properties in synthetic 
models. 

Key words: stars: evolution - stars: fundamental param- 
eters - stars: interiors - stars: late-type - stars: AGB and 
post-AGB - stars: novae, cataclysmic variables 



1. Motivation of the work 

It is well known that the quiescent luminosity of a shell- 
burning star with degenerate core is essentially determined 
by its core mass, without any dependence on the envelope 
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mass, as extensively described in several works carried out 
in the past (e.g. Eggleton 1967; Paczyhski 1970; Tuchman 
et al. 1983). This property is usually referred to as the 
core mass-luminosity (M c — L) relation, although other 
structural parameters may affect the luminosity evolution, 
as indicated in the following. 

Stellar evolutionary calculations actually confirm that 
M c — L relations are generally followed, for instance, by 
(i) low-mass stars during their ascent on the Red Giant 
Branch (RGB) up to the He-flash (e.g. Boothroyd & Sack- 
mann 1988); (ii) low- and intermediate-mass stars dur- 
ing the quiescent inter-pulse periods of their Thermally 
Pulsing Asymptotic Giant Branch (TP-AGB) evolution 
(e.g. Iben & Truran 1978, Wood & Zarro 1981, Forestini 
& Charbonnel 1997), provided that hot-bottom burning 
is not operating (see Blocker & Schonberner 1991); (iii) 
Planetary Nebula nuclei as long as the H-burning shell is 
active (e.g. Vassiliadis & Wood 1994); and (iv) Nova sys- 
tems during their stationary nuclear burning phases (e.g. 
Tuchman & Truran 1998). 

These relations are known to be quite different for 
RGB and TP-AGB stars, and also vary among stars in the 
same evolutionary stage but with different envelope chem- 
ical composition. In view of interpreting these differences 
as the effect of different physical conditions, it would be 
advisable to derive quantitatively the dependence of the 
quiescent luminosity of a giant (RGB and TP-AGB) star 
on basic quantities, namely: core mass M c , core radius R c , 
and chemical composition. 

To this aim, we adopt the same formalism as fully 
described in Refsdal & Weigert (1970), which is based 
on the use of homology relations applied to the case of 
stars with high-density cores surrounded by nuclear burn- 
ing shell(s). The authors demonstrated, for instance, that 
for low values of the core mass (i.e. M c 0.45M©) and 
negligible radiation pressure (i.e. (5 ~ 1) the luminosity 
is expected to depend on the envelope mean molecular 
weight as L oc fi 7 ~ 8 . Such theoretical prediction is found to 
describe extremely well the significant composition depen- 
dence of the luminosity of low-mass stars evolving along 
the RGB, as shown by calculations of full stellar mod- 
els (e.g. Boothroyd & Sackmann 1988; see also Sect. 3.1). 
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However, the results of Refsdal & Weigert (1970) cannot 
be straightforwardly extended to the case of stars evolv- 
ing through the AGB phase, since the condition (3 ~ 1 is 
generally not fulfilled, due to the increasing importance of 
radiation pressure. 

Therefore, our target is to generalize the formalism de- 
veloped by Refsdal & Weigert (1970) for any value of (3 
in the range (0, 1). To do this, we follow the indications 
suggested by Kippenhahn (1981), who first pointed out 
that the quite different M c — L relations for RGB and 
AGB stars are indeed the expression of a unique relation 
modulated by the variation of (3. However, as Kippenhahn 
(1981) explicitly derived the dependence of the luminosity 
only on M c and R c , in this work we will extend the same 
kind of analysis to the composition dependence as a func- 
tion of (3. In this way, we can predict to which extent the 
chemical properties of AGB stars (w ith (3 < 1 ) may affect 
their M c — L relation (see Sects. ^2, and |3.3|) . 

Moreover, given the generality of the method, we can 
also investigate the sensitiveness of the M c — L relation 
to the dominant nuclear energy source operating the H- 
burning shell, i.e. the degree of temperature dependence 
of the relevant reaction rates. This will turn out to be 
significant in view of interpreting the peculiar luminosity 
evolution of RGB stars with initial zero-metallicity, where 
the p-p chain (not the CNO cycle as is usually the case) 
provides most of the stellar energy (see Sect. pTl| ). 

Finally, we invite the reader to refer to the works by 
Refsdal & Weigert (1970) and Kippenhahn (1981) for a 
full understanding of the analytical derivation, since many 
details will be omitted here to avoid redundant repetitions 
and lengthy demonstrations. However, the basic steps will 
be indicated in the next sections. 



The mean molecular weight /x explicitly appears only in 
the equation of state: 



P = Pg 



(3^ 



(1) 



assuming that the total pressure, P, is the sum of the 
contributions of the gas (supposed perfect), Pq, and of 
the radiation, Pr. The quantity (3 is defined as the ratio 
Pq/P. For a fully ionized gas /i = 4/(5A + 3 - Z). 

The parameter eo is related to the rate of energy gen- 
eration by nuclear burning e, which can be conveniently 
approximated as: 



(2) 



with n and v being determined by the rates of the nuclear 
reactions under consideration. Typical values are: (n = 2, 
v ~ 4) for the p-p chain, (n = 2, v ~ 14 -j- 20) for the 
CNO-cycle, (n — 3, v ~ 22) for the triple a-reaction. 

The parameter kq is related to the opacity, which can 
be expressed: 



K P a T b 



(3) 



with the exponents a and b depending on the dominating 
opacity source. Note that in the case k is mostly due to the 
Thomson electron scattering (k = 0.2(1 + X) without any 
dependence on pressure and temperature), we get kq — 
1 + X, and a = b = 0. 

Then, adopting the homology relations presented in 
Refsdal & Weigert (1970) we can express the luminosity 
as a power-law relation: 



LcxM c 5l i? c V 3 Ko 4 £o 5 



(4) 



2. The analytical method based on homology 
relations 

We will describe in the following the analytical procedure 
adopted to derive the dependence of the luminosity on M c , 
i? c , and chemical composition. 

Concerning this latter it can be seen, from the basic 
stellar structure equations, that the effect of the compo- 
sition enters the M c — L relation via three parameters 
(Refsdal & Weigert 1970; Tuchman et al. 1983): 

1. the mean molecular weight /i; 

2. the product of the abundances of the interacting nuclei 
involved in the nuclear burning of hydrogen, eo = X 2 
(for the p-p chain; X is the hydrogen abundance in 
mass fraction) or e = XZqno (for the CNO-cycle; 
Zcno is the total abundance of CNO elements); 

3. a factor k,q expressing the composition dependence of 
the opacity k. 



which, under the assumptions of a fully ionized gas, with 
electron scattering opacity, can be written: 



5X + 3-Z 



(l + X)^(XZ cno ) k (5) 



for dominanting CNO-cycle, or 

■ 63 

L oc M^R s r 2 [ 

c c \5X + 3-Z 



(1 + X) Si (X 2 ) 5s 



(0) 



for dominating p-p chain. 

The exponents (Si, i = 1,5) are the unknown quanti- 
ties to be determined, as a function of (3, of the opacity 
parameters (a, b), and of the energy parameters (n, v). 

The composition dependence is the first being consid- 
ered here, for its particular relevance to the applications 
discussed in the second part of this work (Sect. ||). For the 
sake of completeness, the results for the de pend ences on 
M c and R c are then briefly presented (Sect. 2.3). 
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2.1. The dependence on // 

The dependence on // can be derived treating it as an 
independent parameter to express the homology relations: 

p(r/R c ) oc ^ 
T(r/R c ) oc //* 

P{r/R c ) oc ^ (7) 
L(r/R c ) oc ^ 

Such expressions imply the assumption that ///// = const, 
at each corresponding point (i.e. r/R c — r' / ' R' c ) of models 
with different core radii (i.e. R c and R' c ). Here r denotes 
the radial coordinate of any point inside the region ex- 
tending from the bottom of the burning shell, r — R c , up 
to a point, r = ro, where the variables p, P, and T have 
already significantly decreased, and where L r = L. More- 
over, it is assumed that the other quantities (i.e. p, T, P, 
L, eo, and Kq) do not vary among models at corresponding 
points. 

The exponents 03, (3^, 73, and 63 are the unknown 
parameters to be singled out. To this aim, we need to in- 
tegrate the equations of hydrostatic equilibrium, radiative 
transport, and energy generation, basing on the propor- 
tionality relations given in Eqs. ([?]), and using the expres- 
sions of Eqs. (§) and (|) for the rate of nuclear energy 
generation and opacity, respectively. We then derive: 

P(r/R c ) oc //" 3 
T 4 - b (r/i? c ) oc M «3+ a7 3+5 3 ( 8 ) 
L(r/R c ) oc /i «w+"ft' 

Since there are four unknown quantities, one more rela- 
tion is needed in order to close the system of Eqs. @. 
This is given by the equation of state (Eq. (|l|)), which 
should be expressed in a more suitable form fully explic- 
iting its dependence on (3. For this purpose, we must con- 
sider that (3 is a function of density, temperature, and 
molecular weight. We can then derive the three depen- 
dences as follows. As indicated by Kippenhahn (1981) the 
logarithmic derivatives of (3 with respect to the density 
and temperature are: 

^ ■ ft (?m =-3u-,) ( 9) 



1 



dlnT 



dln P/T,» v 

Similarly, we can make a step ahead and derive also 
the logarithmic derivative of (3 with respect to the mean 
molecular weight: 



dln/3 



= (3-1 



(10) 



Hence, the equation of state can be re-written: 

P oc p- VT 4 - 3/3 (11) 

Note that in the limiting cases (3 = 1 and (3 = 0, we obtain 
the right thermodynamical dependence of the total pres- 
sure when due to the sole contribution of gas and radia- 
tion, respectively. It should also be remarked that Eq. (0) 



is expected to apply for all r/R c . However, since pressure 
changes by several order of magnitudes throughout the 
shell and radiative buffer above, it is far from obvious that 
Eq. (|ll]) does apply throughout the region. In fact it does, 
for two reasons: (1) in RGB stars because (3 ~ 1, and (2) 
in AGB stars because (3 is constant throughout the shell 
and sub-convective layers ("radiative zero" condition). 

At this point, all the necessary information is available. 
Comparing the exponents of Eqs. M), (||), (|TT|), we can 
solve the system of 4 algebraic equations in the unknowns 
«3, #3, 73, and S 3 , yielding: 



13: 



+ 6- 4) 



(v + b- 4)(1 -0) + (1 + a + n)(4 - 3/3) 



j3(l + a + n) 



(12) 



(13) 



{v + b - 4)(1 - /3) + (1 + a + n)(4 - 3/3) 
73 = «3 (14) 



(3[v(l + a + n) - n(v + 6-4)] 
{v + b- 4)(1 - (3) + (1 + a + n)(4 - 3/3) 



(15) 



It is worth noticing that for (3 = 1 we get exactly the 
same results as given by Refsdal & Weigert (1970) (see 
their Eqs. (27) and Sect. II d)), as expected. 

More generally, Eq. (|l5|) allows us to investigate the 
//-dependence of the luminosity in the whole range < 
[3 < 1, so that it is possible to set quantitative predictions 
both for RGB stars (with (3 ~ 1) and AGB stars (with 
[3 < 1). The results are displayed in Fig. |l|, assuming a = 
6 = 0, and characteristic values of the the parameters (n, 
v) corresponding to relevant kinds of energy sources, as 
already indicated in Sect. |^. 

From the inspection of the bottom-left panel of Fig. [l] 
the following features should be noticed: 

— In the case (3=1 and dominant CNO cycle, the expo- 
nent (53 ~ 7, which well reproduces the /i-dependence 
of the luminosity for RGB stars as indicated by cal- 
culations of stellar models (Boothroyd & Sackmann 
1988); 

— In the case (3=1 and dominant p-p chain, we get 
83 ~ 4, i.e. a weaker //-dependence; 

— In any case, for (3 ~ the //-dependence vanishes as 
expected when the gas pressure goes to zero; 

— In the case (3 ~ 0.5-0.8, which are typical values found 
in low-mass AGB stars, and dominating CNO-cycle, 
the //-dependence is not at all negligible. Specifically, 
it results ^3 ~ 1 + 3, a range which is consistent with 
the finding by Boothroyd & Sackmann (1988) (i.e. L oc 
/x 3 ) in their analysis of low-mass (hence with larger (3) 
AGB stars. 

2.2. The dependence on eo and kq 

In an analogous way as described in the previous section, 
we aim now at deriving the dependence of the luminosity 



4 



Marigo P. 



10 
8 
6 
4 
2 




III 

- 


1 1 1 1 1 1 1 


1 1 1 1 1 1 1 

-f 






■I 






a«:f- 






:l ' 






■V '. 






■V -? 






■) cm/ 






r .y y/ - 














Ill 


III 


III 



0.2 0.4 0.6 0.8 1 
P= P g /( p g + p r ) 




'3- 



0.2 0.4 0.6 0.8 1 



P= p g/( p g + p r) 




-2 
-4 
-6 
-8 
10 



1 1 


1 1 1 1 1 1 1 1 1 


III 






-\ \ 






•N _ 






'■.\ 






'•A - 






\\- 






\v 
\\ 




1 1 1 1 1 1 1 1 


1 1 1 1 1 



0.2 0.4 0.6 0.8 1 

P= p g AV p r ) 



0.5 



-0.5 



~i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — r 





i i i i i i i i i i i i i 



0.2 0.4 0.6 0.8 

0= P g /( P 8 + P r) 



Fig. 1. Expected be- 
haviour, as a function 
of (3, of the exponents 
in Eq. (|): Si (related 
to the core mass), 82 
(related to the core 
radius), ^3 (related to 
the mean molecular 
weight), 64 (related 
to the opacity), and 
65 (related to the nu- 
clear reaction rates). 
We adopt a = b = 
for all cases, with 
(n = 2, v = 4) for 
the p-p chain (solid 
line), (n = 2, v = 14) 
for the CNO-cycle 
(dashed line), and 
(n = 3, v = 22) for 
the triple a-reaction 
(dotted line). 



on both eo and kq, already defined in Sect, 
write the homology relations: 



We first 



p(r/R c ) cx Kq 4 



T{r/R c ) cx k p 







Ml 



P(r/R c ) cx K^etf 
L(r/R c ) cx 4*4' 



(16) 



assuming that the functions Ko(r/R c ) and eo(r/R c ) scale 
up by a constant factor at corresponding points of models 
with different core radii. 

Then, integrating the basic stellar equations, express- 
ing the total pressure P as in Eq. (jll]), comparing the 
exponents relative to the four variables p, T, P and L, 
and solving the system of 8 algebraic equations, we finally 
derive the 8 parameters: 

= (4-3/3) 

4 (4-6-i/)(l-/3)-(l + a + n)(4-3/3) 1 ' 
a 5 = 74 = 75 = a 4 (18) 



Si 
S 5 



Ik 



(1 -0) 



(4 - b - u)(l -p)-(l + a + n)(4 - 3/3) 

n(4 - 3/3) + v{\ - 0) 

(4 - b - v){\ - J) - (1 + a + n)(4 - 3/3) 
1 n(4 - 3/3) + u( 1-/3) 

(4 - b - u)(l - /3) - (1 + a + rc)(4 - 3/3) 



(19) 
(20) 
(21) 



Again, for j3 = 1 the above equations yield the same results 
as in Refsdal & Weigert (1970), but now accounting for 
any possible value of (3. 

Concerning the exponents related to the luminosity, 
S4 and £5, their predicted trends are illustrated in the 
bottom-right panel of Fig. [l]. We can notice that: 

— their effect on the luminosity is usually much weaker 
than that produced by except for values of /3 very 
close to zero; 

— in the relevant range of (3 for AGB stars, the predicted 
values for £5 are in a rather good agreement with the 
results by Boothroyd & Sackmann (1988), who quoted 
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L oc Z" n ° 4 basing on their evolutionary calculations of 
AGB models; 

— in the case (3 ~ and a = b = we get 64 ~ — 1 and 
<5 5 ~ 0; 

— both 64 and J5 do not vary significantly adopting dif- 
ferent values of the parameters a, b, v, and n, within 
a reasonable range. 

2.3. The dependences on M c and R c 

In a similar fashion, we can also derive the dependences 
of the luminosity on M c and R c 



L oc M^R S C 2 



(22) 



For the sake of conciseness, we report here only the results 
for 81, and 82 (see also Fig. [j], top-left and top-right panels, 
respectively): 



h = 



and 



na(4 - 3/3) - n(4 - b) - v(l + a/3) 
(4 - 6 - i/)(l - /3) - (1 + a + n)(4 - 3/3) ' 



(23) 



_ /3(3n + 3b + 1/ - 3) - (3 + n)[b + a(4 - 3/3)] 
2 ~ (4-6-i/)(l-/3)-(l + o + n)(4-3j8) ( j 

omitting the formulas for (a,, /3j, 7*; i = 1, 2). The above 
equations for 5i and $2 recover exactly those presented by 
Kippenhahn (1981) setting a = 6 = 0, and those derived 
by Refsdal & Weigert (1970) with /3 = 1. 
In brief, we can outline the following: 

— Si and 82 show opposite trends with (3. However, it 
should be recalled that, though M c and R c are for- 
mally treated as independent parameters, an M c — R c 
relation is expected to exist for highly condensed cores 
which may be assimilated to white dwarf structures 
(Chandrasekhar 1939; Refsdal & Weigert 1970; Tuch- 
man et al. 1983). Considering the homology depen- 
dences of the luminosity only on M c and R c , we get 
the differential equation (see Eq. 15 in Kippenhahn 
1981): 



d In L d In R c 

= 5i + 62 x 



d\nM c 



d In M c 



(25) 



where d In R c / d In M c is the slope of the M c — R c rela- 
tion. In other words, the quantity 82 x dhiR c / d\nM c 
represents the dependence on the core mass via 
R C (M C ). This can be estimated from Eq. (p5|), e.g. 
taking the right-hand side from the results of stellar 
models, and evaluating <5i and 82 from the homology 
relations. 

As for /i, the dependence on M c grows at increasing 
(3. For [3 ~ 1, we get typical relations L oc Mj for 
dominant CNO-cycle, and L oc for dominant p-p 
chain. Again, the former prediction very well agrees 
with the results of evolutionary stellar calculations of 
RGB stars (Boothroyd & Sackmann 1988). The effect 



of a weaker M c -dependence predicted in the latter case 
will be investigated in Sect. 3.1. In both cases (Si ~ 7 
and 82 ~ 4), the luminosity should be just moderately 
affected by changes in the core radius, i.e. the prod- 
uct |<5 2 x d\nR c /d\YiM c \ (calculated with Eq. (p5|)) is 
typically less than unity for /? ~ 1. 
— At decreasing (3 the dependence on M c quickly ap- 
proaches a linear one (8\ = 1 for f3 = 0), whereas the 
dependence on R c tends to vanish (82 = for j3 = 0). 
This result is consistent with the typical flatter slopes 
(~ 1 — 2) of the fitting M c — L relations derived from 
evolutionary calculations of AGB stars. 

3. Some applications 

3.1. The M c - L relation on the RGB: the 
zero-metallicity case 

Evolutionary calculations of low-mass models (see Figs. || 
and |^) indicate that 

1. towards lower metallicities and for given stellar mass 
the luminosity at the tip of the RGB decreases, 
whereas the core mass increases; 

2. at decreasing metallicity the luminosity on the RGB is 
lower for given core mass; and 

3. the slope of the M c — L relation on the RGB is flatter 
for Z = models than for Z ^ models. 

The first point clearly confirms that the luminosity is not 
solely a function of the core mass. To better interpret the 
above trends we can make use of the analytical relations 
obtained in the previous sections. 

Let us consider the M c — L relation presented by 
Boothroyd & Sackmann (1988) 



(11.6M cM ) 7 Z c V 12 



(26) 



which is a fitting formula of evolutionary calculations for 
RGB stars with different metallicities (masses and lumi- 
nosities are in solar units). 

The check the reliability of homology predictions, we 
first compute the exponents Si and £3 with the aid of 
Eqs. d25| ) and (|l5|), and compare them with those given in 
Eq. (|26|). To do this, /3, the opacity parameters, and the 
energy parameters should be specified. We set (3 = 1 - 
which is proved to be a good approximation for all RGB 
models here considered -, and a = b = - for the sake of 
simplicity. The adoption of more proper values for a and 
b would result in very small corrections, as pointed out by 
Refsdal & Weigert (1970). 

The choice of the energy parameters requires some 
comments. It turns out that the bulk of energy produced 
in the H-burning shell is provided by the CNO-cycle for 
both 1M Q RGB models with Z = 0.019 and Z = 0.0004, 
whereas energy production is dominated by the p-p chain 
in the 1M Q RGB model with Z = 0. It follows that typical 
values (n = 2, v = 14 — 16) and (n = 2, v = 4 — 6) should 
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Fig. 2. Left-hand side panel: Evolutionary tracks up to the He-flash of 1 Mq models with different initial metallicities 
(i.e. Z =0.030, 0.019,0.008, 0.004, 0.001, 0.0004 from right to left; taken from Girardi et al. 2000) (thin lines), and 
zero metallicity (thick line; Marigo et al. 2000, in preparation). Right-hand side panel: luminosity and core mass on 
the RGB tip of 1 M Q models as a function of metallicity. 



be adopted in the two cases, yielding Si ~ S 3 ~ 7 — 8 and 
S 1 ~ S 3 ~ 4.0 - 4.7 for RGB models in which the CNO 
cycle and p-p chain dominates, respectively. 

The predictions for the former case (CNO dominated) 
are in excellent agreement with Eq. (26), thus reproducing 
the results of complete stellar calculations of RGB stars 
with dominating CNO cycle. Moreover, it is worth remark- 
ing that points (1) and (2) - mentioned at the beginning 
of this section - are explained as the effect of differences 
in the mean molecular weight, i.e. for given M c and lower 
/x, L is lower. 

The predictions for the latter case (p-p dominated) 
fully explain point (3). In fact, the weaker temperature 
dependence of p-p reactions results in a flatter slope of the 
log M c — log L relation followed by the 1 M Q model with 
initial zero metallicity. A good fit is obtained adopting 



L = (11.6McM) 



1.55 



(27) 



with L and M c in solar units. This power-law relation has 
the same base as in Eq. (|2^), but a different exponent, de- 
rived under the assumption of dominating p-p reactions 
with energy parameters (n = 2, v ~ 5.6). Taking for n the 
value after the first dredge-up as indicated by the calcu- 
lations by Girardi et al. (2000), and letting M c vary over 
the relevant range, we finally obtain the relation shown 
in Fig. |3|, which remarkably well matches the evolutionary 
results at initial zero metallicity. 

To summarise the conclusions, we can notice that, for 
<~ 1, the mean molecular weight determines the luminos- 
ity level of the M c — L relation (i.e. the intercept in loga- 



rithmic plot) for RGB models with similar energy proper- 
ties, whereas the kind of nuclear energy source affects the 
rate of the geometric increment of the luminosity with the 
core mass (i.e. the slope of the log M c — log L relation). The 
latter point explains the fact that the log M c — log L rela- 
tions of the IMq models with Z = 0.019 and Z = 0.0004 
are almost parallel, whereas the relation for Z = runs 
flatter. 

3.2. The composition dependence of M c — L relation on 
the TP-AGB phase 



Let us now apply the relations derived in Sects. 2.1 and 



.2| to TP-AGB stars, in order to estimate the degree of 
dependence of the quiescent luminosity on envelope com- 
position. 

First of all, we need to know how (3 varies during the 
TP-AGB evolution of a stellar model, so that we can di- 
rectly obtain S3 from Eq. (15). To this aim, we make use 



of the definition of radiative gradient: 

i 3 Ijf jPr 



WnacG M r T* 



(28) 



with usual meaning of the quantities. In the context of the 
present study, the radial coordinate r refers to the top of 
the H-burning shell or, equivalently, to some point near 
the bottom of the radiative inert region (extending up to 
the base of the convective envelope). Here L r /M r is nearly 
constant and equal to L/M c , the opacity is dominated 
by Thomson electron scattering, i.e. k = 0.2(A + 1), and 
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Fig. 3. M c - L relation on the RGB. 
Symbols correspond to selected mod- 
els for 1 Mq stars with differ- 
ent metallicities, and mean molecu- 
lar weights after the first dredge-up as 
indicated. Predictions from homology 
relations are shown as solid lines. See 
text for further explanation. 



the radiative gradient, V r , approaches the "radiative zero" 
value of 0.25 (see, for instance, Scalo et al. 1975). Under 
these conditions, and reminding that P r /T 4 oc (1 — 
we can express f3 as a function of the core mass M c , surface 
luminosity L and hydrogen abundance X in the envelope: 

= 1 - 7.956 x 10~ 6 (1 + X) L I L& . (29) 

; M C /M Q 

Figure || illustrates the results of synthetic calculations 
for three TP-AGB models with original solar composition. 
The evolution of the quiescent luminosity (top-left panel) 
is computed according to the prescription presented by 
Wagenhuber & Groenewegen (1998). This gives an accu- 
rate fit to the results of extensive full evolutionary cal- 
culations of the TP-AGB phase (Wagenhuber 1996), ex- 
pressing the quiescent luminosity as a function of the core 
mass, from the first thermal pulse. 

Then, once f3 is evaluated with the aid of Eq. ( p9| ) for 
current values of L and M c (top- left panel) , the composi- 
tion parameters 83 (Eq. |l5|), 64 (Eq. f20| ), and <$ 5 (Eq. plh 
can be calculated (right panels). In other words, we are 
able to predict the current composition dependence of the 
quiescent luminosity as the star evolves on the TP-AGB. 

These simple synthetic calculations are meant to be 
indicative examples, showing the expected sensitiveness of 
the luminosity to possible changes in the surface chemical 



composition during the evolution caused, for instance, by 
convective dredge-up episodes. For the sake of simplicity, 
these events are assumed not to occur in the cases under 
consideration, to avoid the complication due to feed-back 
effects. We will discuss this point below in this section. 

An interesting point to be noticed in Fig. |J is that the 
composition dependence is quite strong during the initial 
(and fainter) part of the TP-AGB phase (e.g. S3 ~ 5 — 6 
in the first stages of the model with the lowest core mass), 
then becoming weaker as the full-amplitude regime is at- 
tained. This trend simply reflects the rate of increase of 
the radiation pressure (i.e. (3 decreases) during the evo- 
lution (bottom-left panel of Fig. ||). Therefore, we would 
expect that, if an efficient dredge- up occurs quite early 
during the TP-AGB evolution, the increase of the mean 
molecular weight could alter the asymptotic approach to- 
wards the full amplitude regime, as otherwise expected for 
unchanged chemical composition. As already suggested by 
Marigo et al. (1999), this prediction, in combination with 
additional effects, may concur to explain the recent re- 
sults by Herwig et al. (1998), who find a steeper increase 
of the luminosity in stellar models with extremely efficient 
dredge-up (with A ~ 1 and larger) already from the first 
thermal pulses. 

Another point to be remarked is that any change in 
the envelope composition is expected to produce a certain 
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feed-back on the luminosity, in addition to the direct effect 
already discussed. In fact, a variation A/x > produces 
AL > 0, since £3 is always positive (we do not consider 
here the extreme case (3 — for which 83 = 0). At the 
same time, we get ST > and 8P < (as > and 
73 < 0), both causing 8(5 < and AL < 0. In other words, 
the effective increase of L due to an increment of /i is 
somewhat reduced with respect to that directly predicted 
by Eq. (filf). Similar effects are produced by variations of 
the other two composition parameters, eo, and kq (see also 
Sect. II g) in Refsdal & Weigert 1970). 

In conclusion, as already suggested by Kippenhahn 
(1981) and confirmed by stellar evolutionary calculations 
(see, for instance Boothroyd & Sackmann 1988), we quan- 
titatively demonstrate that the composition dependence 
of the luminosity of TP-AGB stars is weaker than for 
RGB stars (due to the increasing importance of the ra- 
diation pressure). However, non- negligible effects on the 
luminosity can be driven by significant variations of the 
mean molecular weight, either related to originally differ- 



ent chemical compositions, or caused by the third dredge- 
up. The occurrence of the latter process already since the 
first thermal pulses could significantly affect the luminos- 
ity evolution, in particular, of low-mass TP-AGB stars 
(with large (3). Finally, we remark that the relations pre- 
sented in Sects. |2.l| and 2.2 can be usefully employed to 



improve the analytical description of the luminosity evo- 
lution in synthetic AGB models. 

3.3. The case of zero-metallicity AGB stars 

The last application refers to the TP-AGB phase of stars 
with initial zero metallicity. In Fig. ^| we show the time 
evolution of the surface properties (L and T e g) and of the 
contributions of the nuclear energy sources (He- and H- 
burning shells) for two AGB models of different masses 
(2.5 M Q , and 5.0 M ), taken from Marigo et al. (2000, in 
preparation) . 

As we can see from Fig. |[ the 2.5 M Q experiences 
weak luminosity fluctuations instead of "normal" thermal 
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Fig. 5. Surface and energy proper- 
ties of AGB models with initial zero 
metallicity and masses of 2.5 Mq (left- 
hand side panels), and 5.0 Mq (right- 
hand side panels), since the first ap- 
pearance of the He-shell thermal in- 
stabilities. The fractional contribu- 
tions of shell nuclear burnings to the 
surface luminosity (Ln/L and Ln e /L) 
are displayed, together with that pro- 
vided by p-p chain relative to the in- 
tegrated energy generation from shell 
H-burning (L pp /L H ) 



pulses, whereas the 5 Mq model shows the occurrence 
of rather strong He-shell flashes. Actually, it has already 
been pointed out by Sujimoto et al. (1984; see also Chieffi 
& Tornambe 1984, and Dominguez et al. 1999) that the 
occurrence of thermal pulses in zero-metallicity AGB stars 
is critically dependent on the core mass and abundance of 
CNO elements in the envelope. To this respect, a detailed 
discussion of our zero-metallicity models is given in Marigo 
et al. (2000, in preparation) and will not be repeated here. 

What we simply aim to do in this work is to test 
whether the homology predictions presented in Sect. [2] are 
able to account for the quite different trends in the sur- 
face properties of the two AGB models here considered. In 
fact, as we can notice from the top panels of Fig. [| the 5 
Mq model is climbing the AGB at increasing luminosities 
(and decreasing effective temperatures), whereas the 2.5 
Mq model is clearly evolving downward on its Hayashi 
track. Moreover, in both cases the p-p chain is negligi- 
bly contributing to the nuclear energy generated within 
the H-shell (L pp /£h ^ 10~ 4 ), i.e. the CNO-cycle is the 
dominant energy source. We also report that the third 
dredge-up is never found to occur in these models up to 
the moment at which the calculations were stopped. 

Once (3 is estimated with the aid of Eq. (^) , the expo- 
nents S3, <$4, and ($5 can be calculated assuming a = b = 0, 



and (n — 2, v — 14), the latter being suitable for dominant 
CNO cycle. The luminosity is then derived from Eq. ([5]), 
taking the values of the core mass and composition fac- 
tors (fi, X, and Zqno) from selected models of the 2.5 
Mq and 5 Mq stars (see Figs. || and 0). It turns out that 
in both cases the trend in the luminosity evolution pre- 
dicted by stellar evolutionary calculations is satisfactorily 
reproduced by homology relations. 

In particular, the oscillating behaviour in the surface 
luminosity of the 2.5 Mq model shows up in correspon- 
dence with that of /?, which presents a mirror-like trend. 
The average increase of j3 actually determines the long- 
term decrease of the luminosity, despite of the progressive 
increment of the core mass. Moreover, it is interesting to 
notice that, in the regime of oscillating luminosity, the 
total abundance of CNO elements - at the point of max- 
imum nuclear energy efficiency in the H-burning shell - 
also presents clear fluctuations, reflecting a similar trend 
in the temperature at the same mesh-point. To this re- 
gard, it should be remarked that although for this model 
the CNO abundance in the envelope is zero (as it is not 
changed by any dredge-up episode) , the CNO catalysts in 
the H-burning shell are self-produced, starting from the 
synthesis of primary 12 C via the triple a-reaction, operat- 
ing at typical shell temperatures. 
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Fig. 6. Time behaviour of relevant 
quantities characterising the evolution 
of 2.5 M Q AGB model, namely: the 
gas pressure over total pressure ratio 
(/?), and the mass coordinate (M c ) at 
the top of the H-burning shell; the 
mean molecular weight in the enve- 
lope (fi); the abundances (by number) 
of CNO elements (Zcno) and hydro- 
gen (X), and the temperature (T s h e ii) 
at the point of maximum energy pro- 
duction in the H-burning shell. The 
"true" stellar luminosity predicted by 
full stellar models is compared to that 
obtained from homology relations (ar- 
bitrarily shifted by a constant on the 
logarithmic axis). 



The increase with time of the pre-fiash maximum lu- 
minosity in the 5 M Q model is also well reproduced by 
homology predictions, being essentially determined by the 
rate of increase of the core mass. Finally, we can notice 
that, contrary to the 2.5 Mq model, in this case (3 is de- 
creasing. 

4. Concluding remarks 

A general formalism based on homology relations is pre- 
sented to derive the structural dependences of the quies- 
cent luminosity of shell-burning stars in different evolu- 
tionary phases (RGB and AGB), with different chemical 
compositions, and different nuclear energy sources. 

The reliability of this formalism is tested through sev- 
eral consistency checks. First, we are able to get exactly 
the same formulas as in earlier similar works for specific 
choices of the parameters, e.g. Refsdal & Weigert (1970) 
for (3 = 1. Second, our predictions are found to be in 
good agreement with some basic results of complete stel- 
lar calculations. In particular, as far as the composition 
dependence of the M c — L relation is concerned, we re- 
cover the finding of evolutionary calculations that L oc fi 7 
for RGB stars with Z cno > 0, and L oc /i 3 for AGB stars 
with 0.5M Q <; M c <; 0.7M Q (e.g. Boothroyd & Sackmann 



1988). Moreover, according to our results, the effect on 
the luminosity of TP- AGB stars produced by significant 
composition changes should be larger in the case of low- 
mass stars that experience the third dredge-up after the 
first thermal pulses. 

The case of zero-metallicity giant stars is also investi- 
gated. We show that the particular luminosity evolution 
of RGB stars with Z = can be very well explained by 
considering not only the dependence on the mean molec- 
ular weight, but also the kind of dominant nuclear energy 
source (i.e. the p-p reactions). Moreover, a good reproduc- 
tion of the luminosity trend of AGB models with initial 
zero metallicity is obtained. 

Finally, it is worth remarking that the analytical pre- 
scriptions presented in this work could be usefully em- 
ployed in synthetic evolution models to improve them in 
accuracy, and to test the effects of different physical condi- 
tions, i.e. dominant nuclear reaction rates, opacity source, 
relative contributions of gas/radiation to the total pres- 
sure, and chemical composition. 
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